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Abstract 

A leading twist expansion in terms of bi-local operators is proposed for the 
structure functions of deeply inelastic scattering near the elastic limit x — > 1, 
which is also applicable to a range of other processes. Operators of increasing 
dimensions contribute to logarithmically enhanced terms which are supressed 
by corresponding powers of 1 — x. For the longitudinal structure function, 
in moment (N) space, all the logarithmic contributions of order ln fc N/N are 
shown to be resummable in terms of the anomalous dimension of the leading 
operator in the expansion. 
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a. Introduction: Of special interest in QCD hard processes are "elastic" limits, in which 
the mass of the final state is small compared to the momentum transfer. In such limits we can 
study, in principle, the transition from short- distance to long-distance dynamics. Typically, 
perturbative calculations receive logarithmic enhancements near the elastic limit. Interesting 
and phenomenologically relevant examples include (but are not limited to): the structure 
functions in deeply inelastic scattering (DIS) in the limit x — > 1 where x = Q 2 /2p ■ q, 
with q 2 = —Q 2 the momentum transfer and p the hadron's momentum, the Drell-Yan cross 
section at partonic threshold, z = Q 2 /s —> 1, with Q the pair mass and s the partonic 
center-of-mass energy squared, and the thrust, T in e + e~ annihilation for T — > 1, where 
T ~ (m 2 + m 2 ) / s, with m 2 the masses of the two final-state jets that characterize the elastic 
limit. 

In each of these cases, logarithmic corrections at leading power have been resummed to 
all orders in perturbation theory at leading and next-to- leading logarithmic accuracy Jl] Q ; 
for example, the terms like a™[ln m (l — x)/(l — x)} + , with m < 2n — 1, in the structure 
function F 2 in DIS. Until now, however, there has been no technique to resum higher-order 
logarithmic corrections beyond the leading power of 1 — T, 1 — z or 1 — x. The phenomeno- 
logical importance of these corrections has been highlighted recently in the context of Higgs 
production M. They may also play an important role in renormalon phenomenology [Bj. 

In this paper, we take a step toward the resummation of logarithmic corrections at 
nonleading powers in the elastic limit. We propose that logarithmic enhancements near 
the boundary of phase space can be analyzed systematically through a bi-local operator 
product expansion. Our procedure supplements the light-cone expansion in DIS at leading 
twist, by identifying new operators associated with leading and nonleading behavior of the 
coefficient functions as x — > 1. It is also applicable to processes, such as thrust, for which 
the light-cone expansion is not available. To be specific, we shall consider the example of 
logarithmic enhancements to the longitudinal structure function Fl, but we expect, following 
our comments above, numerous other applications. 

For Fl power singularities at x — 1 are absent. In moment space, the leading power terms 
(in N) are of order a n s ln fc N/N. We shall argue that all such logarithmic contributions can be 
treated systematically within the formalism of large- x factorization and the bi-local operator 
product expansion introduced below. This expansion is distinct from, although consistent 
with, the light-cone expansion, because operators of different dimensionality all contribute to 
the leading power in Q, the difference being only that higher dimensional operators produce 
terms suppressed by higher powers of 1 — x. 

Every inclusive DIS observable is derived as a particular projection of the hadronic tensor 

W, u {p,q) = ^jd 4 y e~** (ptft(O) j v (y)\p) . (1) 

We shall denote such a four-dimensional Fourier transform, with respect to momentum k, as 
FT^ 4) below. For our discussion, it is enough to consider only electromagnetic interactions. 
Generalization to the electroweak case is straightforward. Structure functions are obtained 
via the projections F r = P^ u W^, for r = 2, L, with 

pr = pr = -rf v + \p?. (2) 
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The resummation of logarithmic corrections at leading power in F 2 as x — > 1 may be 
regarded as a consequence of the following factorization theorem || : 

F 2 {x,Q 2 ) = \H 2 {Q 2 )\ 2 J 1 dx' J X X dw J ((a/ - x - w)Q 2 ) V(w) 4>(x') 

= \H 2 {Q 2 )\ 2 J ®V ®<P, (3) 

where (g> denotes the convolution in the longitudinal momentum fraction. This theorem has 
been proven in the axial gauge in Ref. Below we shall define each of the factors that 
enter Eq. (§) in a manifestly gauge independent way H 2 {Q 2 ) is the short distance hard 
scattering function on either side of the final state cut. V is the soft radiation function that 
contains all enhancements coming from on-shell propagation of low frequency partons. It is 
defined as 

V(w) = J dy-e~ iwy ~ pd (0|$t(0,-oo) ^(0,y~v) $ v (y-v,-oa)\0) 

= FT% (0|$t(0,-oo) $e(p,y-v) $ v (y-v,-oo)\0) . (4) 

Here the light-like incoming direction is denoted by v M (+ direction) and the parity reflected 
by v^ 1 (— direction), with v • v — 1. In the second line we have introduced a notation for 
Fourier transforms along the light cone. The Wilson line operator along the v light-cone 
direction is 

& v (x + tv , — 00) = Pexp ^— ig s J ds A^(x + sv)^j . (5) 

<p is the parton distribution function that contains the collinear enhancements from the 
initial state. Its operator definition, for the scattering of quarks, is 

<f>(x')=FT% (p|V(0) jt>$*(0,y-)j(y)\p)®V- 1 , (6) 

where V~ l is defined by V~ l <S> V = 1. The role of V~ x is to remove soft contributions, which 
would otherwise be double-counted in both and V . Finally, J is a function that describes 
the outgoing "current jet", of invariant mass (x' — x — w)Q 2 . The partons in J are almost 
collinear and moving in the light-cone direction. The operator definition of J is 

J ((1 - z)Q 2 ) = FT% p (0|$t(0, -00) V(0) $(y) $ v (y, -00) |0) ® V' 1 . (7) 

The Wilson lines make the functions manifestly gauge invariant. Whenever they do not 
explicitly appear in our expressions henceforth it is to be assumed that they are always 
there. 

b. Operator expansion: Let us now generalize the analysis of the leading contributions 
to F 2 . We shall see that the method not only applies directly to Fl, but also opens the way 
to the analysis of nonleading singular corrections in F 2 , and elsewhere. 

As x — > 1, Fl(x,Q 2 ) depends on three scales: the momentum transfer Q, the hadronic 
mass scale of order Aqcd, and an intermediate scale (1 — xY^ 2 Q, which is essentially the 
mass of the final state. We assume that, although x is close to unity, nevertheless (1 — x) 1 ^ 2 
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^ A-qcd/Q, so that the intermediate scale remains perturbative. At the same time, log- 
arithmic behavior in 1 — x arises from the limit in which the final state becomes massless 
relative to Q. Logarithmic contributions to this limit are associated in perturbation theory 
with the momentum configurations illustrated by Fig. 1. The figure shows a cut diagram 
representation of the scattering process, in which the final state consists of a single "jet" of 
particles, moving in the v^ 1 direction, in the notation of the previous section. Following the 
general analysis of Ref. [7|], in these regions of momentum space the lines of any diagram fall 
into one of four categories. They are either "soft" lines (subdiagram S in the figure), with 
momenta that vanish relative to Q in all four components, "hard" lines, which are off-shell 
by order Q 2 , or "jet-like" lines, with one momentum component of order Q, either in the 
(parallel to the incoming hadron) or (part of the outgoing jet). Hard lines are all con- 
tained in the subdiagrams labelled H (in the amplitude) and (in the complex conjugate) 
in Fig. 1. Note that the hard amplitudes are not specific to any structure function. The 
discussion so far is at the level of the tensor W^ v . 

Because lines in H and are off-shell by the largest momentum scale, Q, we may extend 
the factorization program somewhat, and treat them in terms of an "effective Hamiltonian" , 
as an expansion in products of local fields, 

#: ff = £a(Q^)of(o), (8) 

i 

where \i is the factorization scale and the operators 0\ carry the quantum numbers of the 
electromagnetic current. The functions Cj are the usual coefficient functions that accompany 
the hard-scattering function. For the leading-power analysis of F 2 jlj], we keep only the 
lowest-dimension operator, the original electromagnetic current j u {0) = ip(0) , y u il ! {0), the 
matrix elements of whose product is then factorized into the parton distribution 0, Eq. 
which uses up two of the quark fields in j, and the jet function J, Eq. ([?]), which uses up 
the other two. 

Suppose that P^ v projects out one of the structure functions. Then we can write 

p^w^ = pr ftw ( p \ ^(0)7^(0) HvhMv) \p) 

« FTg(p| ^(0) V(y-i;) \p) ® FT$ xp (0\m $(y) |0> lv . (9) 

Noting that the Wilson lines, Eq. (|), are needed to formulate the factorization in a gauge 
invariant fashion, we see that the leading order factorization of F 2 is thus reproduced. It is 
obvious that this procedure gives a vanishing result when the j's are projected with P£ v in 
Eq. @) to give Fl. All we have to do in this case, however, is to consider the operators of 
the next higher dimension in H e g, Eq. Upon using the equations of motion for both ip 
and ?/>, the following candidates, all of dimension four, can be identified: 

0£ = ^A,<V, 0» d = $D^. (10) 

It is easily seen that and 0<i c have vanishing projection when contracted with P£ v ' . 
They do not contribute to F L but can contribute to F 2 at the subleading level. O^d has a 
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longitudinal projection that is proportional to ip (v- D)ip- After integration by parts this 
term vanishes by the equation of motion of ip at x — 1. We therefore conclude that the only 
operator of dimension four that will contribute to Fl at leading level is 02 a with longitudinal 
projection 

VpO% a = $ p ± i/>. (11) 

Recall that the field ip will be contracted with partons from the incoming state and that it 

is i/i p± that enters the jet function. The reason that the covariant derivative is associated 
with x/j and not on ip comes from dimensional considerations, as explained below. 

To match the above bi- local operators to the general cut diagram in Fig. 1, we need to 
decide which fields inside each 2 will be contracted with partons from the incoming state, 
and which with the outgoing jet. A moment's thought shows that, for an incoming quark 
state, only the quark field if) can be associated with the parton distribution. This is because 
dimensional considerations require C 2 (Q) to behave as Q^ 1 , up to logarithmic corrections. 
The only way to make up this suppression, and derive a result that is not power-suppressed 
in Q, is to integrate over final-state phase space. If fields of dimension 3/2 + d are emitted 
into the final state, the resulting phase space integral behaves as [(1 — x)Q 2 ] d , which precisely 
makes up for the dimensional suppression of C 2 , for d — 1, and indeed for any of the higher- 
dimension coefficient functions in the effective Hamiltonian Eq. No such compensation 
can arise when the extra fields are contracted with lines parallel to the initial state, because 
the initial state is not summed. Consider, then, a general operator Oi in the expansion of 

H e fi, 

Of = Ef V, (12) 

where is some operator of dimension greater than 3/2. The leading power contribution of 
Oi to P^W^ as x — > 1 will be associated with the following factorization, which generalizes 
Eq. (H) directly: 

PrW, u = P? FTf (p\ m 3^(0) Zt/y) ^(y) \p) 

« FTg (p\ m $(y-v) \p) ® FTfl xp Pr(0\ S ft< (0) E^y) |0) . (13) 

The second correlation function is a jet function and we see that at each level in the expansion 
of H c h a new set of jet functions is generated that differ in their contributions by powers of 
(1 — x) according to the dimensional counting rule given above. 

c. The Fl factorization: Next we apply the above analysis to the x — > 1 behavior of Fl- 
Its leading behavior, which comes from O^, is of course leading power in Q, but suppressed 
by a factor of 1 — x compared to F 2 , which is just what we expect. In addition, it is also 
easy to see that for a gluonic state, the lowest-order operator, is even higher order, so that 
the contribution to Fl from gluons is actually finite in the limit x — > 1. The factorization 
formula for the leading corrections of Fl is 

F L (x,Q 2 ) = \H L (Q 2 )\ 2 J'® V®<j). (14) 

The factors V and 4> are the same as in F 2 case, see Eqs. (f|), @. The new jet function J', 
generated by the 2a operator is defined by 
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J'((l - z)Q 2 




(0|$i(0, -oo) p±il>(0) p±ip{y) -oo)|0) ® V~ l . 

(15) 



The conventional normalization factors of Eqs. (|T|), @ have been included. Once the def- 
inition is given in terms of an effective operator, the anomalous dimension of J' can be 
computed in dimensionally regularized perturbation theory with massless partons. Details 
of this calculation will be published in a forthcoming paper. The following points, however, 
are worth noting here: First, one has to be careful to distinguish between the UV renor- 
malization of the operator in J' that corresponds to genuine collinear enhancements, from 
the renormalization of the coupling constant. This is as easy task since the renormalization 
of the coupling is well known. Second, since J' starts at 0(a s ) whereas J, Eq. (0), starts 
at 0{oP s ), the running coupling will mix with the anomalous dimension of J' already at the 
first non-trivial order 0(a 2 ). 

The resummation of the leading corrections to Fl is a consequence of its factorization. 
First, the factorization formula is written in moment space in terms of the Mellin transformed 
factors: 



F L (N,Q 2 ,e) 



( (p ■ Vf ( P • v) 2 2 



1 ~ Tf f Q 2 {P-V? , 2 , 



■ a s (ji 2 



(16) 



with p 11 the jet momentum at the elastic limit. The Mellin transformation of the jet function 
J' is defined without including the overall 1/N, 




(17) 



The factor of 1/N arises because, as seen from the dimensional counting of the previous 
section, the jet J' is of order 1 — x relative to the F2 jet J. It is this additional factor of 
1 — x that shows up as a power suppresion 1/N in moment space. Once this is isolated, the 
renormalization group can be used to resum the logarithmic in iV terms. 

The resummation formula of the logarithmic corrections is derived using the, by now 
standard, methods that are applicable to a large range of processes near their elastic limits. 
Enhancements from soft gluon radiation exponentiate. The evolution of the incoming jet 
can also be written as an exponential. These two factors are common to both F2 and Fl and 
we shall not reproduce them here (for review and details see Ref. 0). We concentrate on the 
differences between the two strucutre functions at leading power in N. The renormalization 
group equation satisfied by the jet function J' is 
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In J' 



Q 2 (p-vf 



Nfi 2 



--7J'K(^ 2 )) 



(18) 



The anomalous dimension jji is novel and characteristic of Fl at the leading power level. 
Up to corrections of 0(a 2 ) it turns out to be M 
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7j'(a s ) = — 

7T 



v_ Cf _ 2Ca _ 4 C(2) (c F _ 9a 



(19) 



The initial condition for Eq. fllq) follows from the lowest order contribution to F^, which 
comes entirely from J'. It is 

J\a s {^)) = C F ^± + 0(a 2 s (/i 2 )) . (20) 

7T 

The final result for Fl may be written in a variety of exponentiated forms depending on 
which anomalous dimensions we wish to include in the exponent and what initial conditions 
we choose for the factors Hl, J', V and 0. All these equivalent forms resum the logarithmic 
corrections at leading power, i.e. all terms of order ln fc N/N . Furthermore, we can exploit the 
fact that apart from their respective final state jets, F 2 and Fl are similar. This means that 
the moments of the coefficient functions in the light-cone expansion C 2 and Cl are connected 
in a precise way, i.e. they differ only in terms that depend on 7j and jj> respectively. In 
dimensionally regularized perturbation theory this relation is 

fl awgyu (21) 

F 2 (N,Q\() 

The ratio of the two structure functions is infrared finite and can be obtained from their 
exponentiated forms. A consequence of Eqs. ( p0|) and ([H]) is that, given C 2 to some order 
in a s , we can readily predict Cl to one order higher in a s . To 0(a 2 ), using the result in 
Eq. (|19"D, we obtain in the MS scheme 

^ ct s C F 
° L ~ ~ AT 

a,\ 2 C F ( Wl , N , (1) 1_ N 



+ n rfL 1 ^\n 2 — -( 1 \ 1 ,'--(3 1 )\n— , (22) 
V vr / 2iV V ^ iVo WJ 2 M ; iV / ' v ; 

with 7# the cusp anomalous dimension Jl| ^[ and A^o = e~ lE . This is in agreement with the 
full 0{a 2 s ) calculation of C L [|,|TJ. 



In conclusion, we emphasize that the operator product expansion presented in this letter 
can be used to provide a systematic treatment of the logarithmic enhancements at the 
edges of phase space at both leading and subleading powers in moment space for a large 
number of hard processes in QCD. Further details and applications will be discussed in 
future publications. 
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FIGURES 



FIG. 1. Momentum configurations that produce logarithmic enhancements near the elastic 
limit x — > 1 in DIS. 
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